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Microorganisms like bacteria can sense concentration of chemo-attractants in its medium very 
accurately. They achieve this through interaction between the receptors on their cell surface and 
the chemo-attractant molecules (like sugar). But the physical processes like diffusion set some limits 
on the accuracy of detection which was discussed by Berg and Purcell in the late seventies. We have a 
re- look at their work in order to assess what insight it may offer towards making efficient, practical 
biosensors. We model the functioning of a typical biosensor as a reaction-diffusion process in a 
confined geometry. Using available data first we characterize the system by estimating the kinetic 
constants for the binding/unbinding reactions between the chemo-attractants and the receptors. 
Then we compute the binding flux for this system which Berg and Purcell had discussed. But unlike 
in microorganisms where the interval between successive measurements determines the efficiency of 
the nutrient searching process, it turns out that biosensors depend on long time properties like signal 
saturation time which we study in detail. We also develop a mean field description of the kinetics 
of the system. 



Berg and Purcell(BP), in their pioneering article [1| 
on "physics of chemoreception" had considered how a 
micro-organism could sense concentration of a chemo- 
attractant molecule (say, X) in its surrounding media. 
They assumed the organism to be a sphere of radius a, 
immersed in an unbounded liquid medium and po be the 
far field concentration of X. A simple example could be 
a bacterium in a dilute sugar (X) solution of local den- 
sity p. The X molecules diffuse and bind to the surface 
of the sphere which is assumed to be a perfect sink for 
X. They solved Diffusion equation, dtp = DV 2 p, in the 
steady state, using spherical co-ordinates centered at the 
sphere. Using the boundary conditions p(r = a) = (i.e., 
fully absorbing surface) and p(r = oo) = po (at far field) 
they obtained the steady state influx of X molecules (J) 
integrated over the spherical surface to be J = AnaDpo. 

We will now briefly introduce a typical biosensor and 
discuss the applicability of the above ideas. A biosen- 
sor is designed to detect traces of specific biochemicals 
present in a carrier medium. It can detect, for example, 
E. coli in drinking water [2[ , hepatitis B surface anti- 
gen present in human serum[3j or pollutants in air [J]. 
The last decade has seen proliferation of such biosen- 
sors [5H7[ in day to day use, mainly due to their, (a) 
quick response time|8], (b) sensitivity to minute amount 
of biomolecules [2] . The particular type of biosensors we 
discuss here are optics based chemical sensors which con- 
verts chemical reactions between GaHIgG (X) and HIgG 
(receptor) molecules into optical signal which is then de- 
tected using fiber-optics technology. In this sensor, an 
optical fiber of radius Ri runs along the axis of a cylin- 
drical chamber of radius R . The fluid containing the 
antigen (X) is injected into the annular space between 
the fiber and the chamber. The surface of the fiber is 
functionalized by putting a certain surface density (cto) 
of antibodies (receptors) on it which serve as the binding 
targets for the antigen molecules. Antigens bind to the 
receptors on the surface of the fiber and absorb evanes- 



cent waves generated by the light carrying fiber. This 
results in loss of intensity carried by the fiber. For our 
purpose here, the absorbance(.A) of the evanescent waves 
[lfj is proportional to the total bound antigen J ad A on 
the fiber surface, where a is the surface density of bound 
antigens. 

Such a system can be described,at the continuum level, 
by reaction diffusion equations [13| |. The X molecules 
bind to the receptors on the fiber surface with a rate uj, 
and surface bound X molecules can also unbind at a rate 
cj„, typically much smaller than the binding rate. The 
values of the kinetic coefficients u b and lo u are unknown 
a priori which we will determine from experimental data. 
The bulk concentration of X is p, the surface concentra- 
tion of receptor-bound X molecules is a and the surface 
concentration of receptors be ctq- Dynamics of p follows 
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S(r - Ri) [p(a - a)ui b - u) u a] (1) 



We will use cylindrical polar coordinate frame where 
p = p(r,<p,z,t) and a = a((p,z,t) with R > r > Ri. 
The second term on the right hand side represents sur- 
face reactions at r = R4. The first term in the square 
bracket describes binding and the second term represents 
unbinding. Dynamics of a follows 
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Here p is the bulk density in the immediate vicinity of the 
surface. These equations can be nondimensionalized. We 
rescale the bulk and the surface densities as p — — , and 



a = — ; the space and time variables as f 
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equations are 
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Superficially the spherical surface of BP is replaced in 
our biosensor by a cylindrical surface but the big differ- 
ence is that our system is confined and the total num- 
ber of antigens is fixed. Thus the steady state here 
corresponds to a state of dynamic equilibrium when the 
binding and unbinding at the fiber surface balance each 
other making both the bulk and surface concentrations 
constant is time. Note that although the surface con- 
centration becomes static the steady state binding (Jj n ) 
and unbinding {J out) fluxes, individually are not zero at 
the surface (see Fi jll2l obtained at different values of uib) 
and at large time Ji n = J ou t = J*. In BP's case, with 
perfectly absorbing surface J, n is given by the surface 
integral of the diffusional current —D J S7pdA onto the 
absorbing surface (of area A). But for our sensor with fi- 
nite binding constant the influx Jj„ = J dAp{a — a)uib is 
computed by integrating the binding term (on the right 
hand side of Eq.2) over the area A of the fiber surface, 
and similarly the outflux J out = J dAuj u a is computed 
from the unbinding term. We consider the initial condi- 
tion where at t = the system is filled up with a fluid 
carrying an uniform concentration of X. At t = the in- 
flux is nonzero as the concentration of X-molecules in the 
vicinity of the fiber surface (p s ) is non-zero. On the other 
hand the outflux is zero at t = because there are no 
bound X-molecules at the beginning. BP has considered 
a perfectly absorbing surface which can theoretically be 
attained in the limit Co — ► o° and w u — > 0. Note that, 
in Eq.l and 2, when <j$ ^$> a the binding term reduces 
to p s aou>b and it appears that we do not need u>b to be 
infinity in addition. But practically <7o is bounded due 
to the finite size of the receptors. BP had approximated 
the receptors to occupy a small area with radius s ~ 10A. 
For our sensor it amounts to about 10 13 receptors cov- 
ering the whole fiber surface. We used this value as the 
maximum coverage a™ ax for Fig HI [2] To compare with 
BP's case we will focus on the binding flux of X only. 

First we adapt BP's general expression [l| for the 
steady flux J = AiraDpo to our cylindrical geometry. BP 
had shown that this flux can be calculated for any shape 
by mapping the steady Diffusion equation V 2 p = to 
the Poisson equation for potential V 2 = 0, in charge free 
space. It can be shown that generally J = AirCDpo where 
C is the capacitance of a conductor with free charge Q on 
its surface. Specifically, C = Q/4>oo, where (p^ is the po- 
tential difference between the conductor and infinity. For 
the sensor the cylindrical fiber is the absorbing surface. 
With a radius Ri = 0.1mm and length L — 50mm (i.e., 
aspect ration 500) it is as good as an one dimensional line. 
For a line charge density A, extending from x — —L/2 
to L/2 the expression for the potential <b{z) along the 
perpendicular bisector, z distance away from the center 

of the line charge, is (biz) = fc — In —, , 



k = I/47reo is the Coulomb force constant. Using the 
approximation R4/L «lwc get, <p(Ri) ~ ^ ln(2L/i?,) 
and C ~ %kx(2L/Ri). Thus J B p ^ 2'wDLp hx{2L/R i ). 
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FIG. 1: Semi-log plot of binding (Ji n ) and unbinding (Jout) 
fluxes (both scaled with the steady state value J») as a func- 
tion of time, at fixed D and po(~ O.OOlmg/ml). Values of cot 
and w„ are same as those found through Fig[3](to be discussed 
later). Ji„ and J out are plotted separately because their scales 
of variations are very different, unlike that in Fig[2] 



To compute the binding flux J;„ we have to numeri- 
cally time evolve the dynamical equations (Eq.3,4). First, 
to get realistic values for the kinetic coefficients ujb and 
lo u we use experimental data on surface adsorption a(t) 
versus time from Ref[9(, obtained at two widely differ- 
ent initial bulk densities a) pn — O.OOIma/mi and b) 
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FIG. 2: Semi-log plot of binding (Ji n ) and unbinding (J ou t) 
fluxes (both scaled with the steady state value J*) as a func- 
tion of time. Here we used the same po, D and ui u as in Fig[l] 
while Wj, was increased 100 times to reach a steady state value 
J t comparable to Jbp ■ We explain later why this comparison 
is not very useful due to difference in the boundary conditions 
in the two problems. 



where 



po = Q.lmg/ml, and possibly at different surface density 
of receptors. Note that the non dimensionalized equa- 
tions Eq |3l4l do not explicitly scale with antigen (X) den- 
sity po and therefore these data sets can be treated as 
independent. Despite the wide difference in po the satu- 
ration times (to) in the two cases were similar (the sym- 
bols in FigfS]). This could be rationalized by noting that, 
in case-b the fiber was soaked in the receptor solution for 
two hours while for case-a it was soaked for a very long 
time (about 16 hours). From this information we inferred 



that in case-a a% = a"k 



while for case-b Cn < cjn 



We choose D = 10 5 cm 2 /sec typical of diffusion of small 



molecules in water [111, |12[ (BP also took the same D for 
their estimates). We had to determine ujb,u u and a^ by 
matching our numerical results (from Eq.3,4) with the 
temporal profiles of ait) and the ratio <J% a j<J h oa - We con- 
verged to uj h = 0.75 x 10~ 5 pm 3 J sec, uj u = 0.35 x 10~ 2 / sec 
and (Tq = 0.014pg/mm 2 . These numbers for ui u ,u>b ap- 
pear reasonable when compared to the reaction-diffusion 
processes on bacterial membrane (13J |. 

As mentioned earlier we used cylindrical polar coordi- 
nate system to discretize the space. Uniform binning was 
used along z and 4>; while r coordinate was binned non 
uniformly such that the volume of each bin (rdrdcpdz) re- 
mains constant. Reflecting boundary condition was used 
at the walls of the cylindrical chamber, by ensuring zero 
currents at the boundaries. We used an uniform distri- 
bution of X molecules in the bulk as our initial condition, 
i.e., p(t = 0) = po and a(t = 0) = 0. 
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FIG. 3: Surface density of bound antigens a(t) versus time, 
y-axis is scaled with the saturation value Uoo = a(t — > oo) 
since we do not know the proportionality constant connect- 
ing experimentally measured absorbance A and a. The sym- 
bols represent experimental data. In (a) po — O.OOlmg/ml 
and in (b) po = 0.1mg/ml. The solid lines are from our 
numerical integration of Eq.3,4. For (a) we chose ao = 
0.08pg/mm 2 i.e., the maximum possible surface coverage. 
It turns out that a reasonably good match with the two 
experimental a(t) profiles and with a^/a^ were obtained 
for oj b = 0.75 x 10" 5 /im 3 /sec, lj u — 0.35 x 10" 2 /sec and 
Oq = O.OUpg/rnrn 2 , i.e., about 5.5 times less than the maxi- 
mal coverage. 



We then compute the flux J,- n (t) = J dAp(ao — cr)u>b, 
which is the binding term in the right hand side of Eq.2, 
integrated over the cylindrical fiber surface, as a func- 
tion of time. Interestingly, Ji n (t) goes through a minima 
before it saturates to J* (see Fig. 1,2). We will explain 



the origin of this non-monotonic behavior later when we 
study the dynamics in detail. In Fig.l the steady state 
flux J* is much lower than Jbp, while in Fig. 2 it is com- 
parable. But J* and Jbp depends on different set of 
parameter values. Both of them are steady state prop- 
erties, but J* depends on u!b,uj u ,ao and po while Jbp 
depends on D and po . This difference arise from the dif- 
ference in the boundary conditions of a confined versus 
an unbounded system. Therefore the comparison is not 
fare. J* can be calculated by setting the left hand side 
of Eq.4 to zero and using mass conservation, which will 
be discussed later. 

So far we had implicitly assumed that the microorgan- 
ism can sense the ambient po by measuring the influx (J) 
of X molecules. But BP had also considered the realistic 
possibility that they can infer po by measuring the state 
of occupation of its surface receptors i.e., density of re- 
ceptors that are bound to X molecules. In fact this is 
the recipe which most practical biosensors employ. For 
example, in our particular sensor a{t) decides the inten- 
sity of optical adsorption. In BP's theory a bacteria can 
sense its a(t) in response to local po and decide to move 
towards or away from the chemo-attractant or the chemo- 
repellent, respectively. But for a static biosensor a(t) 
can only increase towards a saturation. Since a system 
takes some time to attain saturation, this measurement 
process is inherently slow compared to the measurement 
of instantaneous flux. On the other hand measurement 
of any instantaneous variable is prone to fluctuation er- 
ror where as long time observables like ait — > oo) are 
more dependable. So the challenge is either to reduce 
the saturation (waiting time) time or choose an optimum 
time interval T over which an instantaneous variable like 
J(t) or a{t) should be measured (so that AJ/J or Aa/a 
is small). BP had correctly concluded that a bacteria 
must employ the second strategy since it has to rapidly 
change its direction of motion based on comparison be- 
tween its successive measurement of cr(t). BP had es- 
timated T ~ lsec for E. coli bacteria. Recent findings 
14j show that bacteria has a very efficient mechanism 
for amplifying the minute signal generated by binding of 
external sugar molecules to its receptors. It has the capa- 
bility of detecting 0.1% percent change in the attractant 
density and that too over four orders of magnitude of 
sugar concentrations. Ref[15[ has shown, that for a par- 
ticular type of biosensor flux detection could be a supe- 
rior method compared to measuring long time saturation 
properties. For our sensor, we now investigate in detail 
how saturation time of the sensor varies in response to 
po and how it can be steered by choosing <j$. 

First we will discuss a simple Mean Field (MF) limit 
of the dynamics. In the MF approximation we consider 
the surface concentration to be uniform over the surface 
of the fiber and the volume concentration to be uniform 
through out the bulk. Let Vq be the volume of the an- 
nular space and Aq be the surface area of the fiber. At 



t = 0, N = p V~o an d later N = p M Vo + cr M A , where 
Pm = PkA^) i s the mean field density (denoted by sub- 
script m) of X molecules and a M — a M (t) is the corre- 
sponding surface density of the bound X molecules. In 
the nondimensional form we have 



where p A 



p M = 1 - aa M , 
and a 
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bulk density p- M being homogeneous and slaved by a M 
(via EqjS]) we need to consider only the equation of mo- 
tion for the surface reaction, namely Eq|3J Substituting 
for p M , from Eqj5l into Eq|4j and simplifying, we get 
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where Ai = ckDj,, A 2 = [(1 + a)u}b + Co u ] and A3 = 0)5. 
Integrating this Equation we get 
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where A/? = \J\\ — 4A1A3. Inverting the above equation, 
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For any set of parameter values, it can be shown that Xr 
is always real and we also have X2/X11 > 1- As a result 
tanh - (X2/X11) is always a complex number. Using the 
standard property tanh - x — coth - x — «7r/2, we can 
rewrite Eq|9]as 
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This formula gives excellent fit to the numerical data (not 
shown here), obtained by integration of Eq.3 and 4, at a 
high value of the diffusion constant, D = 10~ 3 cm 2 /sec. 

The steady state solution (a s u ) can be obtained either 
by setting a M =0 in Eq.6 or from the r — > 00 limit of 
Eq|U We getV = ail [A2 - X H ]. 

The mean field approximation will fail if diffusion is 
not sufficiently fast compared to the time scale at which 
surface binding reactions cause a depletion in the anti- 
gen concentration (p). In such a scenario the spatial in- 
homogeneity in p (along r) takes a long time, compa- 
rable to the saturation time of the sensor, to homoge- 
nize. A better understanding can be gained by compar- 
ing the time scales of the three processes: diffusion (to), 
binding (£5) and unbinding (t u ). We get the individual 
time scales from Eq.l, by comparing each term on the 
right hand side with the left hand side. For example, 
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FIG. 4: Radial density profile p as a function of r — Ri in 
the bulk, at three different times: right after start (square), 
at saturation (circle) and at some intermediate time (trian- 
gle), (a) and (b) differ in the parameter N s /N (which is 
proportional to 00 at fixed po)- The values of D, w&, ui u were 
obtained through Fig[3] Transition from mean field to non- 
mean field type density profile occurs as we go from (a) to 
(b) by increasing N s /N . But note that, at a fixed po, the 
fraction of antigens (X) remaining in the bulk can be reduced 
(consequently the bound proportion can be increased) by in- 
creasing N 3 /N . This is desirable for making the sensor more 
sensitive, specially when po is small. 



p ~ DV 2 p gives, by dimensional analysis, t D ~ j^-. 
Similarly, t^ 1 ~ ^ and t" 1 - ^. Here we have 
assumed R = R — Ri to be the only relevant length 
scale. For diffusion, this is the spatial scale of density 
inhomogeneity. Now, tb and t u are the time scales over 
which density inhomogeneity are created near the fiber 
due to the surface reactions, while £d is the time inter- 
val during which such inhomogcncities are ironed out. 
Therefore, mean field approximation requires diffusion 
to be a faster process, i.e., to <C tb,t u . These inequali- 
ties yield the criteria ^^ < 1 and ^^f^ < 1. The 
first inequality suggests that mean field approximation 
will be correct at high D or low a values. We have 
verified these conditions numerically by looking for den- 
sity inhomogeneity p(r) during the transients, in the nu- 
merical solution of Eq.3, 4 (see FigQJ. For example, For 
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FIG. 5: Semi-log plot of binding (Jj„) and unbinding (J ou t) 
fluxes as a function of time in the mean field regime (i.e., high 
po and low oo with N s /No — 0.05). The symbols (circles and 
squares) are obtained from numerical integration of Eq.3,4 
while the solid lines are the corresponding mean field results. 
Here po = 0.1mg/rnl,ao = 0.01figm/mm 2 and the values of 
cjf, and ui u are the same as those in Fig.l (i.e., the values 
estimated from experimental data). 
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FIG. 6: Semi- log plot of saturation time to (sec) versus 

the density of receptors ao(pg/mm 2 ). The inset shows sat- 
urated signal (Too versus <jq. We fix po at a very low value 
O.OOlmg/ml to test the sensitivity of the sensor. The values 
of W;, and u, are those estimated before. Both t and (Too sat- 
urate at high values of oo, much beyond the maximum surface 
coverage (0.08pg/mm 2 ) considered here. We have explored 
the seemingly unrealistic <7o > a™ ax regime here because it 
may be possible to increase <to by choosing smaller receptor 
molecules in another system. 



N s /No = 1,D = 10~ 5 cm 2 /sec the density remain uni- 
form through out, at all times. But when oq is increased 
by choosing N s /Nq = 10, strongly inhomogeneous p(r) 
appears (i.e., MF theory fails). Now in addition if D is 
hiked p{r) becomes homogeneous again (graph not shown 
here). The second inequality suggests, along with high 
D and low &q, we also need high po- Then only both 
tb,t u ^> to can be satisfied. We have verified this con- 
dition on po along with similar conditions on Wf, and lo u 
resulting from the inequalities. FigjS] shows a compari- 
son between numerical solution of Eq.3,4 and mean field 
results for p = O.lmg/ml and <jq = 0. 01 \xgmj 'mm 2 , i.e., 
at high po and low ero . At these parameter values the in- 
flux Ji n does not go through any minima The reason why 
the influx Ji n goes through a minimum in Fig |ll2l is now 
clear from Fig2b, which shows when p(r) becomes inho- 
mogeneous (in the non MF case) the p(r) in the vicinity 
of the fiber undergoes a dip (triangles) before it becomes 
uniform (circles) at late times. In the mean field regime 
the minima is absent because the p{r) in the vicinity of 
the fiber decreases monotonically in time as is clear from 
Figlli. 

We now study the general case when MF theory is in- 
valid and thus we have to depend on numerical integra- 
tion of Eq.3,4. Our numerical curves for a versus time, 
shown in FigOJ could be fit to exponential functions like 
c(£) = 0-00(1 — exp _t / T °), allowing us to estimate a satu- 
ration time scale tq and the saturated value a^ (plotted 
in Fig |6l7l) . To and <Joq depend on both oo and po- 

The aim of Fig |fj] and [7] is to identify the regimes where 
saturation time To can be reduced and saturated signal 
(Too can be maximized. For FigEl Pa has been held fixed 



300 



250- 

o 
H 

200 h 



150 



- 


1 


- 


- 





- 




• 






0.08 




• 




II 


• 






.y 


• 




. 


0.06 


• 




• 




O 


• 




• 


- 


0.M 


/ 




• 
• 




0.02 






• 


" 




..•• 




• ~ 


- 





» • • •* - ^ 




• 




0.001 0.01 0.1 p 1 




- 


1 


1- 



0.001 



0.01 



0.1 



FIG. 7: Semi- log plot of saturation time to (sec) versus the 
density of X molecules po(mg/ml). Here cto is held fixed at the 
maximum surface coverage (0.08pg/mm 2 ) in order to maxi- 
mize the signal. The inset shows that even at this maximum 
surface coverage the saturated signal, a x drops drastically at 
low po. 



at a low value, O.OOlmg/ml, while for FigEl oo is fixed 
at a™ 11 * = 0.08fj,g/mm 2 . The insets of both the figures 
show that signal can be enhanced either by increasing po 
or (To, which result in decrease or increase of To, respec- 
tively. Of course at high p a strong saturated signal can 
be achieved within a short saturation time, but the sen- 
sitivity of a sensor is tested when p is small which we 
will focus on below. For low po , uq should be maximum 



to maximize the signal, even at the cost of higher wait- 
ing time. Operating near maximum receptor coverage is 
also necessary as Fig |5] shows that the nonlinear response 
starts to increases near this point. But for moderate and 
high po, we should choose moderate cr such that To is 
not so high and the signal is strong enough. This may 
appear analogous to the conclusion of BP where with just 
a fraction of the cell area (~ 1/1000) covered with recep- 
tors the steady flux could be as high as Jbp/2, where 
Jbp is the maximum flux with the fully absorbing sur- 
face. But the assumption behind this derivation was that 
the inter-receptor distance is much much greater than the 
receptor size. In Fig|6l at 1/10-th of the maximum sur- 
face coverage (i.e., at er ~ 0-0lp,g/mm 2 ) the signal a^ is 
much weaker compared to that at a™ ax . This again high- 
lights the difference between our confined system and the 
steady state behavior of Berg and Purcell's unbounded 
system. 

In summary, we examined the applicability of Berg 
and Purcell's ideas to real sensors. In general it turns 
out that a flux based sensor is more efficient than one 
which depends on long time signal. The flux in our 
sensor also shows unexpected time variation which 
results from competition among different time scales and 
the extended nature of our system. Another interesting 
observation is that even at realistic diffusion constant, 
mean field theory works when po is high and Co is small. 
In general, nonspecific binding of X molecules on the 
fiber surface can cause complications but for the system 
we have chosen here nonspecific binding was verified to 
be negligible. Further, the surface reactions need not be 
first order, which we have assumed here. We checked 



that consideration of second order binding kinetics does 
not give any new exotic behavior (eg, oscillations etc) 
but changes the quantitative values of saturation time. 
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